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Abstract 

We obtain the Boltzmann equations for leptogenesis including decay and scat- 
tering processes with two and three body initial or final states. We present an 
explicit computation of the CP violating scattering asymmetries. We analyze 
their possible impact in leptogenesis, and we discuss the validity of their ap- 
proximate expressions in terms of the decay asymmetry. In scenarios in which 
the initial heavy neutrino density vanishes, the inclusion of CP asymmetries in 
scatterings can enforce a cancellation between the lepton asymmetry generated 
at early times and the asymmetry produced at later times. We argue that a 
sizeable amount of washout is crucial for spoiling this cancellation, and we show 
that in the regimes in which the washouts are particularly weak, the inclusion 
of CP violation in scatterings yields a reduction in the final value of the lep- 
ton asymmetry. In the strong washout regimes the inclusion of CP violation in 
scatterings still leads to a significant enhancement of the lepton asymmetry at 
high temperatures; however, due to the independence from the early conditions 
that is characteristic of these regimes, the final value of the lepton asymmetry 
remains approximately unchanged. 



1 Introduction 

Leptogenesis [1] provides a very attractive scenario to generate the baryon asym- 
metry of the Universe. Indeed, the existence of heavy right handed neutrinos 
is strongly motivated by the see-saw mechanism [2] introduced to generate the 
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observed light neutrino masses. The Majorana nature of the right handed neu- 
trinos imphes violation of lepton number, the large values of their masses imply 
that the out of equilibrium requirement can hold when they decay, and moreover 
these decays can be CP violating due to the phases that are generally present 
in the Yukawa couplings with the lepton doublets. Since the three Sakharov 
conditions [5] can be met, the question then is a quantitative one, i.e. whether 
the asymmetry generated in these scenarios is large enough to account for the 
observed value. 

In recent years, quantitative analyses of leptogenesis have become more and 
more sophisticated, taking into account many subtle but significant ingredients, 
such as several washout processes [H El [71 |8], proper subtraction of on-shell 
heavy neutrino intermediate states [9], thermal corrections to particle masses, 
couplings and decay asymmetries [9], spectator processes [10l[Tl], flavor effects 
[lailllllllllllllBlinillHlllS], and the possible effects of the heaviest right handed 
Majorana neutrinos A^2,3 |121 [20l [2T1 [22] (for reviews of the most recent results 
see [23]). 

The aim of this work is to discuss in some detail the inclusion of CP viola- 
tion in scatterings, and of processes involving two and three body initial or final 
states. In Section 2 we derive the Boltzmann Equations (BE) involving these 
terms and obtain a general parameterization for them. In Section 3 we present 
a detailed computation of the CP violating asymmetries in scatterings in zero 
temperature field theory. These asymmetries have been considered before in the 
limit in which they are proportional to the CP asymmetry in decays [71 I16j . 
and we discuss here the validity of this approximation. In Section 4 we stress 
how in the thermal leptogenesis scenario, when one starts with a vanishing initial 
abundance for the right handed neutrinos and assumes that no other interactions 
besides the Yukawa couplings can produce them, washout processes are essential 
to generate a lepton asymmetry. This is due to the fact that the lepton asym- 
metry produced at early times (mainly through scatterings mediated by off-shell 
A^i and processes producing real A''i's) tends to be compensated by the oppo- 
site sign asymmetry produced at late times (when A''i disappearance processes 
become the dominant source). The inclusion in the BE of CP violating asym- 
metries in scatterings is important to make this cancellation complete. In the 
presence of washouts, the cancellation can be partially avoided because washout 
processes erase slightly more efficiently the asymmetry produced at early times 
than the asymmetry produced at later times. However, when washouts are par- 
ticularly weak the cancellation remains effective, and in this case the inclusion 
of CP violation in scatterings leads to a sizeable suppression of the final lepton 
asymmetry. In the strong washout regimes washout processes attain thermal 
equilibrium, thus erasing any dependence from the earlier conditions and in par- 
ticular from the large asymmetries generated by the scattering processes. Since 
at lower temperatures the contributions from scatterings are highly suppressed 
and the lepton asymmetry is essentially generated through decays, for the strong 
washout regimes the inclusion of scattering CP asymmetries leaves the final re- 
sults approximately unchanged. 
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2 The Boltzmann equations 



We consider thermal leptogenesis scenarios with hierarchical heavy neutrinos Ni, 
i.e. with masses M2^3 ^ Mi. We will focus on the evolution of the lightest one, 
that in this section will be denoted simply by N, and we will ignore the possible 
contributions from the heavier neutrinos A''2,3 |22| except for their virtual effects 
in the CP violating asymmetries. 

Let us first introduce some notation. We denote the thermally averaged rate 
for an initial state A to go into the final state B (summed over initial and final 
spin and gauge degrees of freedom) as: 

7^ = 7(^-i?), (1) 
and the CP difference between the processes for particles and antiparticles as 

A7^ = 7^-7|. (2) 

Particle densities are written in terms of the entropy density s, i.e. Ya = Ua/s 
where is the number density for the particle a. To simplify the expressions 
we rescale the densities Ya by the equilibrium density Y^*^ of the corresponding 
particle, defining = Ya/Y^'^, while the asymmetries of the rescaled densities 
are denotecfl by Ay a = Va - Ua- 

The difference between a process and its time reversed, weighted by the den- 
sities of the initial state particles, is denoted as 

n m 

[^-i?]=(nyaj7^-(ny^.)75, (3) 
i=i j=i 

where the state A contains the particles ai,...,a„ while the state B contains 
the particles bi,...,bm- We will consider only processes in which at most one 
intermediate state heavy neutrino N can go on the mass shell, and in these cases 
a primed notation 7^ will refer to the rates with the resonant intermediate state 
(RIS) subtracted (and similarly for [A <-> B]'). Accordingly, a process A ^ B 
of this kind will be divided into a RIS subtracted (off-shell) piece, and a second 
part corresponding to on-shell N exchange: 

lB=l B+1 B- (4) 

In the simple case when only 2 <-> 2 scatterings are considered, the on-shell part 
is just 

7"^^ ^ 7i^i?f , (5) 

where B^ denotes the branching ratio for N decays into the final state B. (As 
is discussed in section 12. 2| the inclusion of 2 ^ 3 scatterings, and in general 
of processes of higher order in the couplings, implies that eq. ([5]) needs to be 
generalized.) 

^The notation adopted here differs from the notation used in [111114] in which the symbol ya, rather 
than A?/a, was used to denote the asymmetries of the rescaled densities. 
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In the following two sections we will write down the BE for the evolution 
of the density of the heavy Majorana neutrinos N and of the asymmetry for a 
generic lepton flavor i. We will first consider in section 12.11 the leading terms 
involving the Yukawa couplings (generically denoted as A) of the neutrino N 
to the Higgs boson H and to a light lepton doublet £j. In section [2.21 we will 
include the additional contributions arising from processes involving also the top 
Yukawa coupling ht and the gauge interactions. Accordingly, the contributions 
to the evolution equation for the N density will be split in two parts: 

Yn={Yn)^ + {Yn)^^, (6) 

where we have introduced the notation Y = szHdY/dz, with z = M^/T and 
H{z) being the Hubble rate at temperature T. In eq. ([6]) {Yn^ ^ includes the 

contributions of terms up to O(A^), while {Yn^ includes the contributions up 

to 0{\^h'l) or 0{\^g^), with g a generic gauge coupling constant. 

The contributions to the evolution equation for the density of the lepton flavor 
i will also be split in different parts: 



yr^^ = {yL)^+{n)^^+{Y,) (7) 



sphal 

with Yl^ = 2Y£. + Yg- , where the factor 2 comes from summing over the densities 
of the two gauge degrees of freedom in ii , and the inclusion of the density Y^. for 
the right-handed lepton is required because the L-conserving charged lepton 
Yukawas can transfer part of the asymmetry to the right handed degrees of 
freedom. In eq. d?]) (Yl-^ ^ includes the contributions of terms up to O(A^), 

while (Yl^^ includes the contributions up to 0{X'^h^) or 0{\'^g'^). 

The last term ( Yl - ) represents the change in the lepton densities due to 

V ' I sphal 

electroweak sphalerons, which are the only source of baryon number violation. 
Although the precise rates of sphaleron effects are hard to estimate, one knows 
that they are efficient below T ~ 10^^ GeV and that they leave unchanged B — L. 
Moreover, sphalerons generate the same change in the baryon number of each 
generation, in fact one has that {Y^L^sphai = [Y^b) sphal /"i, where Yal, = ^L, - 
Yi^ and Iab is the baryon asymmetry to entropy ratio. Hence, it is convenient 
to write directly an equation for the quantities Y^i = Y^b /3 — Y^Li which do 
not depend on the sphaleron rates. By subtracting from eq. ([7]) the analogous 
equation for Y^. and by subtracting again the result from the equation that 
describes the evolution of the baryon asymmetry, (Ya_b)/3 = (5^A_B)sp/ia«/3, one 
obtains 

Ya^ = -{Yal,)^-{Yal,)^^. (8) 
2.1 1^2 and 2^2 processes 

Let us briefly sketch the way in which the BE are obtained, by considering flrst 
the lepton (flavor) number violating 2^2 scatterings mediated by N exchange. 
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that is £iH ^ ijH (with j ^ i) and iiH <^ ^jH. The BE describe the evolution of 
the particle densities at any given time, but since the Universe is expanding, the 
temperature is decreasing and the particle densities are changing, only processes 
that to a sufficiently good approximation can be described by (effective) contact 
interactions can be consistently included. However, the 2^2 scatterings we are 
considering can proceed through the exchange in the s-channel of an on-shell A^, 
and in this case they are characterized by a time scale that can be comparable 
to the expansion rate of the Universe (at T ~ Mi). Therefore they cannot be 
approximated by a contact interaction, and must be treated with care. The usual 
way to deal with this is to separate the off-shell part of the scatterings from the 
on-shell piece, i.e. 

e^H ^ ,l_,H osljH .gS 
'e^H ~ I ijH ^ I IjH' ^ > 

and similarly for iiH ^ ^jH. Then, at each given time, the evolution of the 
density of the lepton flavor Lj is determined by the following (instantaneous) 
reactions: 

i) N ^ iiH decays occurring at a rate proportional to the density Yat of the 
right handed neutrinos, and iiH ^ N inverse decays; 

ii) off-shell 2 2 scatterings j'l^^ and j'^^^ involving only the exchange of 
virtual A^'s (in the first process N can be exchanged in both the s and t 
channels, while in the second process only in the s channel); 

Hi) 2 <-> 2 scatterings j£^£^ and 7^|^ with N exchanged in the t and u-channels 
(in these cases no RIS can appear, and hence there are no on-shell contri- 
butions to be subtracted). 

The corresponding BE then reads: 

sub / . \Nt 



where 



yl,),^, = w^m, (11) 

sub 
2^2 

Nt 
2^2 



^^^02^2 = J2i^jH^iiH]' + Y.i^,H^iiH]', (12) 
{^L.)Z. = T.{[fiH^iii,] + {l + 6,,)[HH^idj]}. (13) 



It is important to remark that while eq. (jlip contains processes of 0{X'^), both 
eqs. (|12p and (jl3p contain only non-resonant scatterings, that are O(A^). How- 
ever, while in a first approximation the contributions in eq. (1131) may be neglected, 
the inclusion of the off-shell contributions of eq. (jl2p is mandatory. This is be- 
cause the CP asymmetries of the subtracted rates are of the same order than 
the CP asymmetries of decays and inverse decays, and therefore neglecting them 
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would yield inconsistent results. We can now subtract from eq. (jlOp the analo- 
gous equation for Yi. and write separately the source and washout contributions 
to Yal, = Yl, - Iz/as: 

(iAL,)^=(iAL,)%(iAL,)J. (14) 

At the leading order, the source term receives contributions from the 1^2 
decays and inverse decays in eq. (jlip and from the off-shell parts of the 2 <-> 
2 scatterings in eq. (jl2p . while the CP asymmetries of the t- and tt-channel 
processes in eq. (fT3]) . being of higher order in the couplings, can be neglected. 
We can then write the source term as: 



/. \s /. \s /. \s sub 

(lA.O, = (^A..),^, + (^A.O,^,- (15) 



where 



V i jj^t I 

In order to proceed we now use an important relation that states that the CP 
asymmetries in the off-shell scatterings are, to leading order in the couplings, 
equal in magnitude and opposite in sign with respect to the CP asymmetries of 
the corresponding on-shell scatterings: 

Ay^ ^ -A7-i (18) 

To derive this relation, we first substitute the definition of the off-shell rates in 
eq. (HI) to get l^ij^ = - A7°*;g and then we use the fact that the CP 
asymmetry of any process is always of higher order in the couplings with respect 

to the corresponding tree level process [23]. Thus the CP asymmetries A'^ln 
and A7^^^ of the full 2 <-> 2 scatterings are ©(A^). On the other hand, the CP 
asymmetries of the on-shell pieces are of O(A^). This can be seen by writing 
them in terms of eq. ([5]) to obtain 

A7 ^ = A7^ + 7Ar - ^lA ■ (19) 

"ftot 

In eq. jtot is the total decay rate, and we have used A7^ = A7^ = 
— A7^ where the first equality follows from the CPT relation j{A —>■ B) = 
7(5 A), and we have approximated at leading order B^ ~ B^. This shows 
that A7°^;g (and hence A7';g) is of the same order in the couplings as the CP 
asymmetry in decays A7^ (that is 0{X^)). Therefore, up to 0{\^) corrections, 
the contributions of the off-shell rates in eq. (jl7p can be written as 

(lA..) - -2A7;r^ E + Bf^s) . (20) 
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Here B^^jj represents the branching ratio for the decay N — > ijH. At the order 
in the couphngs we are working here J2j i^^fjH + ^fu) ~ ^ therefore the 
r.h.s. of eq. (|2Up can be further simphfied to —2A"f^^. After summing up the 
two contributions (fT6|) and ([20|) . the source term (Yl^^ ^ in eq. ([15]) becomes 
proportional to — 1. This is in agreement with the general condition that 
no asymmetry can be generated in thermal equilibrium, and constitutes a check 
that, at this order, all the relevant contributions to the source term of the BE 
have been included. ^ 

The washout term (VALi^ ^ in eq. (|14p contains the terms proportional to 
the light particle asymmetries, and is the sum of three different contributions 
obtained by subtracting from eqs. pT|) - (fT3|) the corresponding equations for Y^. . 
It can thus be written as 

(\w / . \w / . \ «),sub / . \ w,Nt 

Yal.)^ = {Yau)^^^ + (lA..),^, + (^AL.),^, . (21) 

After linearizing in the CP asymmetries and in the asymmetries of the normalized 
densities, these contributions read: 

(^^^OLa = + AyH)7^'', (22) 

^^Ofr--E[(A.. + A,.)(y|S + 4-) 

+(Ay,^.+Ay^)(7'|5-7l;H)]> (23) 



{yAuY^^^^ = - E [(1 + ^^^^)(^y^. + A2/£, + 2Ay^,)7|| 



(24) 

One can estimate the off-shell parts in eq. (|23p by introducing a subtracted 
propagator for N, or alternatively eliminating the 7' by means eq. dH). 

Finally, to compute the lepton asymmetry we also need to solve for the evo- 
lution of the heavy neutrino density y^v that appears in the source term eq. (|16p . 
To the leading order in the couplings, the corresponding BE reads 

= E {[^^-^ ^ ^] + [^^-^ ^N]}^ -{y^ - 1)7^^2^ (25) 
j 

where 7^^^ = J^jiltjH + ^l^i?) thermally averaged two body N decay 

rate, and we have approximated the leptons and Higgs particle densities with 
their equilibrium values. 



2.2 1^3 and 2^3 processes 

We are now ready to generalize the above procedure to include processes involving 
the Higgs Yukawa coupling ht to the right handed top quark, which is sizeable. 
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Processes involving gauge bosons can be included in an entirely similar way and 
will be mentioned later. 

The inclusion of AL = 1 processes like N ^ iiQt decays and inverse decays 
(where Q is the left-handed quark doublet and t the right-handed top singlet) 
and scatterings mediated by Higgs exchange like Nii ^ Qi, follows along lines 
analogous to those presented in the previous section. For the evolution of the 
heavy neutrino density we obtain 

(iiv)^^ = -iVN - 1) [7^^' + > (26) 

where 

7r^-E(<Q.+7gg,-) (27) 
j 

is the contribution from decays into three body final states, and the contribution 
from Higgs mediated scatterings is 

2^2 _ V- ( Ne, Ni, NQ NQ Nt , Nt_\ .38^ 

where the first two terms in the sum correspond to the contributions from s- 
channel Higgs exchange, while the other four terms (that at leading order are all 
equal) correspond to the contribution from t and u channel Higgs exchange. 

Regarding the evolution of the lepton asymmetries, the derivation of the BE 
is more delicate because, besides the inclusion of the CP violating asymmetries 
in 1 <-> 3 decays like ^ iiQt and in 2 ^ 2 scatterings like Nii ^ Qt, the 
asymmetries of various off-shell 2 <-> 3 scatterings, that contribute to the source 
term at the same order in the couplings, should also be included. Accordingly, 
the term (y^.^ in eq. ([7j) can be written as 



with 



11. 1^3 = [N^ iiQt] + [Qt ^ Nii] + [Nt ^ Qi,] + [NQ ^ tii] ; (30) 



2^2 
sub 

2^3 



\ sub ^ — , , _ , , _ _ , 

Yl^^^ = ^ Q^^'^ + ^ tii] + [ijHt ^ Qii]' 

+ [ijQt w iiH]' + [ijQ ^ iiHt]' + [ijt ^ i,HQ]'} 



+ J2 { ^ W]' + [^jHQ ^ tii]' + [ijHi ^ Qii\ + [IjQt ^ iiH\ 
j _ 

+ [ijQ ^ iiHt]' + [iji ^ iiHQ]' + [Qt ^ iiHij]' + [Qt ^ iiMj]'} 
+ Y.('^ + 6,j)[Qt^iiHij]'; (31) 



, Aft 



{Yl,)Z^ = J2{i^jQ^^ W + [ijH ^ iiQt] + [ijHt ^ iiQ] + [ijHQ ^ iit]] 



8 



+ J2{[QtH ^ ij^ + [QtH ^ Ijii] + [QH ^ (jiit] + [tS ^ IjiiQ]} 
j 

+ ^(1 + 6ij) { [iS ^ £iQ£j] + [QtH ^ litj] + [QH ^ liUj]} . (32) 



As in the previous section, the asymmetries of the off-shell 2^3 rates in eq. (j3ip 
can be estimated by relating them to the asymmetries of the corresponding on- 
shell parts. However, for 2 <-> 3 scatterings the definition of the on-shell part is 
more subtle, because after a real is produced in a collision, it has a certain 
probability to scatter before decaying. Let us consider for example the process 
ijHQ —5- tii. The contribution to this process from the exchange of an on-shell 
N corresponds to the production process ijH N followed by the scattering 
N + Q ^ tii mediated by a Higgs in the t channel. Processes of this kind can 
generally be written as AX — > Y where A denotes a possible state to which N 
can decay. The corresponding on-shell rate then is 

where we have introduced the quantity Py^ that is the probability that the 
heavy neutrino N will scatter with X to produce Y. Processes in which the 
on-shell can disappear only by decaying (as for example ijH iiQt or ijt 
iiHQ) can generally be written asA^BorasX^ BY where both A and B 
denote possible final states for N decays. The corresponding on-shell rates are 

7°^B = InPb, (34) 
7°^^y = 1§yPI (35) 

Note that because of the fact that in the dense plasma can suffer inelastic scat- 
terings before decaying as is described by eq. ([33]) . the quantities in eqs. 
and (I35p differ from the usual notion of branching ratios at zero temperature. In 
particular, scattering rates should also be included in normalizing properly the 
decay probabilities, and the quantities then will denote the general probabil- 
ities that the heavy neutrino A'" contained in state a ends up producing state b. 
In the case under discussion, we have for example 

pN _ lijH r)N_ _ 'iiQt 

" 7a// ' " 7a// ' 

pNii ^ 'Qt pNQ ^ hlj pNt ^ ^Q^^ j-ggN 

7a// ' 7a// ' lall ' 

with similar definitions for the probabilities of the CP conjugate processes. The 
probabilities are normalized in terms of the sum of all the rates, that reads 

ill.H + ll,H+li^Qt + ^m+lQt+^Qt +^U+^U+^1,Q + H^q)- (37) 

i 

To the order in the Yukawa couplings that we are considering, the unitarity 
condition for the sum of the branching ratios of A^ into all possible final states 
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J2y = 1 is then generalized to J2x v Py'^ ~ ^- That is, the probabiUties for 
all the possible ways through which N can disappear add up to unity. 

To include the new sources of CP asymmetries, we now need to subtract from 
eqs. (p0|) - (j32|) the analogous equation for Yi_. We obtain 



(iALjf/ = (iALJt„3+(^AL,) 



2^2 



s, sub 



2*^3 



(38) 



where we have neglected the CP asymmetries of the 2^3 processes with N 
exchanged in the f-channel since they are of higher order in the couplings. For 
the first term in the r.h.s. of eq. (|38p we have 



ALJU, = iVN + 1) [AY;^q, + AY;t + ^^NQ _ 
2^2 



Qt 



(39) 



After eliminating the subtracted rates by writing their CP asymmetries as minus 
the CP asymmetries of the on-shell rates and keeping terms up to ©(A^/i^), we 
obtain for the second term in eq. (j38p 



s, sub 
2^3 



-2 



Qt 



(40) 



Note that at order ©(A^/i^) eq. ^ reads 

, s sub 



(^AL,) 



2*-*2 



(41) 



where the sum of the probabilities Yl,j \Pe^H + Pi Rj unity. However, the 

first term in eq. (j40p conspires with eq. (j4ip to yield the correct behavior for the 
source term involving /S.^f_^, that turns out again to be proportional to {un — 1). 
By summing up eqs. p^ . ([39]) . (|l0]) and (jH]) we obtain the final expression for 
the source term that holds at 0{X'^hf): 



(^alJ 



s 

I+II 



N 



I) 



Nt 



(42) 



Regarding the washouts, the contributions from eqs. (|30|) - (|32|) after subtract- 
ing the analogous equations for can be written as 



\w / . \w / . \w, sub / . \ 



ui, Nt 
2^3 



where 



(^AL,) 



1^3 
2^2 



{AyQ - Ayt - Ay(,^)'y^^Q^ + (Ayg - Ayt - yNAyi^)i 



Qt 



+ {Ayq - ywAyt - Ay^j7^* + {ywAyQ - Ayt ~ Ayi^)-f_^ 



NQ 



iQl 



(43) 



(44) 
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+ {Aye, - Ayg + Ay, - Ayn - Ay,,) ( + 7^, + 7'&, 



+ \{-^yij + Ayg - Ayt - Ay^^ - Ay^) x 



+ i^yej - ^yq + ^yt - ^yn - Ay^.)7^^|^^. 

+ {Aye^ + Ayq - Ayt + Ayn - ^yiihfs,. 



+ ( Aye, - Ayg + Ay, - Ay^ - Ay4)(7§| + 7§f + 71""^ 



+ E [(^y'^j + - + Ayff - Ay£,)7^ 



QtH 

• .f. 



+{^yej + AyQ + Ayt - AyH - Ay£.)7 



tH_ 



(45) 



+ 



+ (Ay^. - AyQ + Ay^ - Ay^^ - Ay^J (7^^ + 7||) 
E {(1 + Sij)i-Aye^ + Ayg - Ay, - Ayn - Ay,J (jj^^^ + 7gf + 7g2J } • (46) 



We can now summarize (and slightly generalize) the procedure for writing 
the BE for leptogenesis. Ignoring processes with more than one N, one can write 
them as 



iAr = -(yiv-l)E7r, 



A,B 



{YAL,r = {yN - 1) E [^^(^) - ^^(^)] ^^r, 



(47) 
(48) 



(lAL,)" = E [^^(^) - ^^(^)] E ^vai - y7 E ^y^i (49) 

where in the first two equations the states A and B contain only light standard 
model particles, while in the last one for the case in which no on-shell intermediate 
N is allowed in the process A ^ B, one has simply 7'^ = 7^. In these equations 
Li{A) denote the i lepton flavor number of state A while n^i = or 1 counts 
the number of A^'s contained in state A (and similarly for Li{B) and ub)- To 
avoid double counting in both equations the sums are restricted to Li{B) > 
and Li{B) > Li{A), where the first restriction avoids double counting the CP 
conjugate processes with li in the final states, and the second restriction avoids 
double counting the time reversed processes where the number of ii in the initial 
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state is larger than in the final state. In the equation (j49p for the washouts 
the Oj's and 6j's denote all the particle species with non- vanishing asymmetries 
contained respectively in states A and B. Let us note that in the cases we 
have been discussing all the CP violating asymmetries appearing in eq. (08]) have 
Li{B)—Li{A) = 1, and thus the expression for the source term can be accordingly 
simplified. However, this is not true for the washout term since some processes 
with Li{B) — Li{A) = 2 also contribute. 

Several contributions that we have included for consistency and for com- 
pleteness, for practical purposes can be neglected without affecting sizably the 
numerical results. This is the case for example for all the non-resonant 2 <-> 3 
washout terms. Also the inclusion of N decays into three body final states has 
been carried out mainly for the sake of completeness, that is to account for all 
the processes of the same order in the couplings, and also to incorporate consis- 
tently those 2 <-> 3 scatterings for which the on-shell piece involves precisely a 
1 — > 3 decay, as for example IjH — > iiQt. However, the overall impact on 
quantitative results of the three-body decay CP asymmetry Aj^q^ (as well as 
the contribution to the washouts of 3 — > 1 inverse decays) are rather small. This 
is because while e.g. A^^q^ involves the same CP violating phase as A'yf,^, it 
also has a significant suppression factor arising from three body phase space. The 
explicit expression for the three body decay rate is presented in the Appendix. 

Following the same procedure outlined above, it is possible to include in the 
BE other relevant processes, such as those involving the gauge bosons I?, 9|. 
With all the subdominant terms neglected and with the effects of the gauge 
bosons included, the simplified expression of the BE for the evolution of Ia? 
reads: 



where A = Wi or B for SU(2) and U(l) bosons respectively, and a sum over all the 
gauge boson degrees of freedom is understood. In eq. (f50]) (as well as in eq. 
below) we have neglected three-body decays like '^^ui- involving the gauge bosons 
because, similarly to the three-body decays involving the top quarks, they are 
suppressed by phase space factors and give negligible contributions, and we have 
also neglected the contributions to the washouts from 2 <-> 3 processes. The 
simplified expression for the evolution equation for the charge Ia^ = ^a_b/3 — 
^ALi that is conserved by sphaleron interactions, is: 



Yn 



[yN-l) Id +ltop +1A 



(50) 



where the term involving the gauge bosons is defined as 




(51) 




(52) 



with 




(53) 
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+ E 



(1 + 5,,)iAyi^ + Ay,^. + 2AyHH^'^ + (Ay,, - Ay,j7ji^ 



(54) 



+(y7vAy,, - Ayg + Ayt)j^^^ - {yN^Vq - ^Vt - ^Vt.hu? 
-{-VN^yt + At/q - Ayi.)-fQl - {yN^ye, + Aynh^f^ 

HyN^VH + ^yt,hAh + (^^^ + ^y^J^^v 

Before concluding this section one more remark is in order. To solve the set 
of coupled BE for the different lepton flavors it is necessary to express all the 
normalized particle asymmetries Ay^ and Ayn in terms of the relevant charge 
asymmetries Ya^ (and use the relation Ayg — Ayj = Ayn/'^)- In doing so one 
has to take into account the constraints on chemical potentials enforced by the 
reactions that, in the particular temperature range in which leptogenesis occurs, 
are faster than the Universe expansion [11^ [25] . Although we have written the 
equations for a generic asymmetry Ya^ of one of the three flavors, it is important 
to keep in mind that when no lepton Yukawa couplings are in equilibrium (T > 
10^^ GeV), one has just an effective one-flavor equation for the asymmetry in 
the lepton i coupled to the lightest heavy neutrino N . When just the r Yukawa 
coupling is in equilibrium (10^ GeV < T < 10^^ GeV) two flavor components are 
relevant, £r and £a, where £a is the component of the lepton doublet produced 
in N decays which is orthogonal to ir [12]. When also the /u Yukawa coupling is 
in equilibrium (T < 10^ GeV) the three flavor components of £ are completely 
projected out by the Yukawa interactions, and the full set of equations for the 
three YJ\^, with i = e, fi, r, is needed. 



3 The CP asymmetry in scattering processes 

We will now study the CP asymmetries in 2 — >■ 2 scattering processes involving 
the top quark and with a Higgs exchanged in the s or in the t channels. These 
asymmetries were considered previously in [7] and [16] where some arguments 
were given in support of an approximate equality between the scatterings and 
the decay asymmetries, as for example: 



Qt 



(55) 



In this section we will explore the validity of this kind of approximations. 

The Lagrangian for the Yukawa interactions relevant for the computation of 
the associated CP asymmetries, written in the mass eigenstate basis of the heavy 
neutrinos Na{a = 1, 2, 3), of the charged leptons {i = e, r) and of the quarks, 
reads: 

Cy = -ha Na£i H'^ -htQtH + h.c. . (56) 
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(a) 





H 





Figure 1: Diagrams contributing to the CP asymmetry in Qt N^ii scatterings. 



Here H = {H^, H^)^ is the Higgs field, with H = iT2H*. Like in the usual case 
of the CP asymmetries in decays, the CP asymmetries in scattering processes 
arise from the interference between the tree level and one loop amplitudes. They 
can be computed by explicit evaluation of the corresponding loop integrals, or 
just by using the Cutkowski rules that give directly the absorptive part of the 
Feynman diagrams. 

The vertex contribution to the CP asymmetry in the scattering Qi <-> N^ii 
with the Higgs exchanged in the s-channel arises from the interference between 
the diagrams in fig. [l^a) and Wib). Expressed in terms of the CP difference 
between the squared invariant amplitudes \A4\'^ — |7V4p this contribution reads 



\MiQi 



6 Pl • P2 P'l ■ P2 I 




Ml M^ + Ml 
-Ml 



.(57) 



Here 6 is the step function {9{x) = 1 for x > and otherwise), s is the squared 
center of mass energy, pi,P2, p'l and p'2 are the momenta of t, Q, ii and respec- 
tively, and Mg is the mass of Nj^ (note that for (3 = a there is no contribution 
since the product of the relevant Yukawa couplings is real) . The overall factor of 
6 corresponds to the summation over the gauge degrees of freedom. Note that 
the term proportional to 6{s — M|) appears due to the possibility of performing 
a new cut in the one-loop graph in fig. W^b) involving the A'^ and lepton lines. 
Since can go on-shell only when the center of mass energy is sufficiently large 
{s > AfJ) this contribution is relevant only for temperatures not much smaller 
than M^. Our results hold in the zero temperature limit; in particular we take 
all the particles, except the heavy Majorana neutrinos, to be massless. At high 
temperatures, finite temperature effects can induce non negligible corrections to 
our expressions. In particular, when Mi < Mh{T) + M£.{T) (implying that de- 
cays and inverse decays are blocked [9j) we expect that thermal masses will also 
have the effect of suppressing the scattering CP violating asymmetries. 

Due to crossing symmetry, the CP asymmetry for the processes N^i Qii 
and NaQ t£i in which the Higgs is exchanged in the t-channel can be obtained 
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from the previous result by replacing the Mandelstam variable s hy t. Note that 
for massless quarks t < so that 0{t — Mj) = and hence in this case no new 
cut is present. 

The cross sections are obtained by integrating the modulus squared of the 
invariant amplitudes: 



a 



64^2(^^+^2)2 7 |pi 



iPlI 



\M\^dn[, (58) 



where pi = (£'j,pi), p'^ = {El,p[) (i = 1,2) are the momenta of the initial and 
final particles respectively, and the cross section for the CP conjugate process a 
is defined in the same way. In terms of the cross sections, the CP asymmetry of 
the s-channel scattering is: 

[aiQt ^ NJ,) - a] (vertex) - - — 



^Im[A»c,A*^(AU)^„ 
9(3 -M. 



in (At A), 




Unlike what happens in the case of decays, the asymmetry a — ci of the scattering 
rates now depends on s, so that the convolution necessary to obtain the asym- 
metry A7^^ of the thermally averaged rates doesn't lead to a simple analytical 
expression, and has to be performed numerically. 

For the CP asymmetry in the i-channel scattering there is no simple analytical 
expression. Also, the usual infrared divergence appears in the limit t — > 0, which 
can be regularized by replacing the factor \/t coming from the (massless) Higgs 
propagator by l/{t — rnj^), and using here the Higgs thermal mass. 

The CP asymmetry in scatterings coming from the wave function piece (in- 
terference between the diagrams in figs. U^a) andlU^c)) turns out to be always 
the same as the CP asymmetry for the decays: 



^^NleSwave) _ Aj^^fj{wave) 



eiiwave), (60) 



where the decay asymmetry is defined in the usual way as ej, = A^^_'^/jj^° 

The ratio between the CP violating scattering asymmetries and their approx- 
imate expressions in terms of the asymmetry in decays derived in ref. [16j, are 
shown in figure [2] as a function of t/mE The results for processes with the 
Higgs exchanged in the s and in the i-channels are presented separately. Two 
illustrative values for the ratio of the two lightest Majorana neutrino masses have 
been used: M2/M1 = 5 and M2/M1 = 100. The effects of have been ignored 



^For definiteness, we neglected in this plot the contribution to the wave part proportional to 
Ma{\'^ \)af)i keeping the one proportional to Af^(AtA)^ct, in which case vertex and wave contributions 
become proportional to the same combination of couplings, see ref. [26) . 
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Figure 2: Ratio between tlie CP violating asymmetry in scatterings A'j]^^_/2j^g, and the 
corresponding quantity for decays Aj^fj/l^^^j, vs. T/Mi for two ratios between the masses 
of the two hghtest heavy neutrinos M2/M1 = 5 and M2/M1 = 100. The contributions from 
s and t-channel Higgs exchange processes are shown separately. 



(this would correspond either to the case M3 ^ M2 or to the situation in which 
the complex phase in the combination of Yukawas associated to A^3 is particularly 
suppressed). It is apparent that this ratio starts to deviate from unity already 
for T > M2/IO and that deviations of few tens of percent can appear for T ap- 
proaching M2. Since the relevant temperature for leptogenesis in these scenarios 
is typically 0.1 < T/Mi < 10, the approximations adopted in ref. [3 [IB] should 
be good if M2/M1 ^ 10, while some corrections appear for milder hierarchies 
at temperatures T > M2/IO. It is also easy to show analytically, starting from 
eq. (|59p . that a factorization for the vertex part, analogous to that of the wave 
part, is obtained in the limit of Mi and T <C M2. 

Regarding the scattering processes with gauge bosons, such as N£ HA, 
NH £A or — > iH, the associated CP asymmetries can be obtained in a 
similar way, computing the interferences between tree and one loop scattering 
amplitudes. One significant difference is that now box diagrams are present, in 
which the gauge boson is attached to a lepton or Higgs in the loop of the vertex 
like diagrams, leading to more complicated expressions. The absorptive parts 
can be obtained using Cutkowski's rule and new cuts appear, but again only 
for s > M|, so that in the hierarchical limit the factorized expression holds for 
T <^ M2, as was the case in the scatterings with quarks discussed before. 

The fractional contributions to the source terms in the BE for Yj\. arising 
from decays and scatterings are shown in fig. [3] using the factorized expressions in 
terms of the decay asymmetries. We present separately the results for two body 
decays Fd = Aj^^^j/ A^, s-channel Fh^ = A7^^/X]A7, t-channel Higgs 
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exchange processes Fh^ = 2A'j^^/J2^7 ^^'^ the gauge boson contributiorH 
Fa = {A^i^J + A70 + E At where 

^ A7 ^ At^Ih + A7,% + 2 A7,7 + A7iy^ + A7,Tf + A7^| , (61) 

while the contributions from decays into three body final states are always neg- 
ligible and are not shown (see the appendix). These fractions are independent 
of the value of the neutrino Yukawa couplings adopted and, in the hierarchical 
limit in which the factorization is valid, they are also independent of the value of 
M2/M1. From fig. [3] it is seen that scattering CP asymmetries are the dominant 
source term for T > 2Mi, and hence can play an important role in the early 
leptogenesis phase|l| 
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Figure 3: Fractional contribution to the source terms from the CP asymmetries in decays (D) 
and in scatterings (s and t channel Higgs exchange and scatterings involving gauge bosons) 
for M2 S> Ml, obtained in the zero temperature approximation. 



4 Results 

In fig. m we display the results of the integration of the BE adopting M\ = 
10^^ GeV and mi = v^{\^X)xilMi = 0.06 eV (where v is the Higgs VEV). 
To avoid complications with r-flavor effects that are active in this temperature 
regime, we also assume a flavor 'aligned' situation in which the lepton doublet 
ii to which A^i decays has no r flavor component, that is Kr = |(^iKt)P = 

■^The contribution of gauge boson scatterings has been estimated using the expressions given in [9]. 

^When thermal masses are taken into account, at very high temperatures (T > 7 Mi) the condition 
Mh{T) > Ml + M(.{T) is met and the decay H — > Niii can occur. Since the asymmetry for this 
decay has a large enhancement from thermal effects [SI , in this temperature regime actually the Higgs 
decays would become the dominant source of the lepton asymmetry. 
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(and similarly for ii) . The dashed line corresponds to the case in which the con- 
tributions of the CP scattering asymmetries to the source term are ignored (see 
ref. [llj for more details), while the solid line depict the results obtained by in- 
cluding the CP asymmetries of the scattering processes, adopting their factorized 
expression (using the exact expressions would slightly modify the asymmetries 
for T > Ml, but the final values would be almost unchanged). The dotted line 
is the asymmetry that would result had the initial density of N be the equilib- 
rium one. It is apparent that for T > Mi the scattering processes have a large 
effect in the production of the lepton asymmetry. This example corresponds to 
a case of strong washout, with rhi 3> 10~^ eV. In particular, we see that in this 
case the washouts affect the evolution of the lepton asymmetries up to 2; ~ 10. 
Hence, even if at early times large additional sources of CP violation are present, 
late washouts turn out to be decisive in determining the final asymmetry, which 
ends up being equal to the one that would be obtained had one started with 
Y/v = This also means that in this regime the final asymmetry becomes 

essentially independent from the conditions at early times. 
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Figure 4: Evolution of the B — L asymmetry as a function of z = Mi/T with the CP 
asymmetries in scatterings neglected (dashed line) and included (solid lines). The scattering 
asymmetries are approximated by using the asymmetry in decays. This figure corresponds 
to a strong washout regime with rhi = 6 X 10-2 eV. For comparison, the value that would 
be obtained starting with an equilibrium N density is also displayed (dotted line). 

We recall here that a crucial point in thermal leptogenesis is that if the 
washout processes were switched off completely during the whole leptogenesis 
phase (and if the dependence of the CP asymmetries on the temperature that 
is induced by thermal effects was also ignored) the inclusion in the BE of the 
sources of CP violation from scatterings would yield a zero final asymmetry [16] . 
This can be seen by writing the different ALj = 1 scattering CP asymmetries in 
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terms of their approximate expressions, i.e. 



A7f ^ ^ ^A7,T^. (62) 

Then, the source term in the BE for Y^- can be rewritten just in terms of the 
decay asymmetry ^'yf.H as 

{Y^.)s^-{yN-l)^ E (63) 

^^^H AL,=1 



By using the relation 



(64) 



where the sum in the denominator in the r.h.s. is over all the processes 7y ^ for 
which |ALj| = 1 and over all the flavors j, and defining the flavor asymmetry 
£\ = A'yf_fj/ J^jiltjH + ^i^^) finally rewrite the source term as 



{YAjs^-{yN-l)e\J2j^''. (65) 

X,Y 

Combining now the above expression with the BE for Y^r we obtain 

^N-—r- — (66) 

In the absence of washouts the r.h.s. of the equation above would then vanish, 
and in the approximation in which e\ is taken as independent of the temperature, 
the quantity Yfq—Yi^-/e\ would hence be constant. Then, for thermal leptogenesis 
scenarios in which the density and the lepton asymmetries vanish initially, this 
quantity will just be zero, showing that the Y^. asymmetries generated at early 
times are erased at later times as the A''i disappear by decays or scatterings. 
Then, as was pointed out in [9], any effect that breaks this cancellation, as 
for example a dependence of the CP asymmetries on the temperature, could 
be numerically important. The cancellation will no longer hold also if some 
washout processes are particularly efficient at temperatures T > Tq, where Tq 
is the temperature at which the lepton asymmetry changes sign (note that this 
type of washouts could yield an enhancement in the final asymmetry, while in 
general late washouts at T < Tq always tend to reduce its final value). Of course, 
in the cases when the heavy states are produced through other processes not 
related to the ones giving rise to the CP asymmetries, such as via scatterings 
involving heavy right-handed W or additional Z' bosons, or are produced non- 
thermally via e.g. inflaton decays, the leptogenesis initial condition y/v(0) ^ 
would directly prevent the cancellation of the lepton asymmetry. 

In other words, the origin of the cancellation can be understood as follows. 
At any time there are three kinds of possible sources for the lepton asymmetries: 
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off-shell scatterings, processes producing real A^s and processes in which A^s are 
destroyed. In general, the lepton asymmetry produced by the off-shell scatterings 
is twice as large, and with opposite sign, as that associated to processes in which 
real A'"s are produced, but that associated to processes in which real A^'s are 
destroyed depends on the N density Yj^. If Y/v equals the equilibrium density 
Y^^, the asymmetry of A'' production and destruction processes are equal, and 
hence the sum of the three asymmetries cancel. On the other hand, if Yat < 
the asymmetry produced by off-shell processes dominates, while the opposite 
happens if Y^ > Y^'. This is why the sources of lepton asymmetry in the BE 
are just proportional to Y^v — ■ Now, when the factorization of the scattering 
asymmetries in terms of the decay asymmetry holds, one finds that the sources of 
the lepton asymmetries are just ei times the sources of the A^ density (assuming 
that no other processes besides the Yukawa couplings produce A^s). This means 
that if ei is constant (T independent) and we ignore the washout processes, the 
total integrated change in the lepton asymmetry will be just ei times the total 
change in Y^r. Hence, if this last vanishes, as is the case when the initial condition 
is that of vanishing A^ density, the final leptonic density would also vanish. 
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Figure 5: Same as Figure S] but in a regime of weak washout, with m,\ = 10 ^ eV. 



Clearly, the impact of including CP scattering asymmetries in the BE is qual- 
itatively different in the strong and in the weak washout regimes. In the weak 
washout regimes (corresponding to values of m-i < 10^'^ eV) the effect of late 
washouts is negligible, the asymmetry is strongly affected by the cancellation and 
thus its final value turns out to be rather sensitive to the amount of washouts in 
the early phases. This is illustrated in the example in fig. O that corresponds to 
the value rh\ = 10"^ eV. In this fi gure we compare the evolution of the asym- 
metries in the two cases when the CP asymmetries of the processes involving 
the top- Yukawa and the gauge interactions are included or are left out (all scat- 
terings are in any case included as sources for A^ production). It is apparent 
that when the sources of CP violation from scatterings are included the effects 
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of the cancellation strongly reduce the final asymmetry obtained. Since only 
weak washouts are present, the cancellation remains quite effective and the final 
value of the asymmetry is rather small. In this regime, the final asymmetries 
are also much smaller than the asymmetries that would result starting with an 
equilibrium density for the A^'s since in that case no asymmetry is generated at 
early times, and no cancellation can occur. 
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Figure 6: Same as Figured] but in a regime of intermediate washout, with rhi = 10 



Finally, fig. [6] gives an example with an intermediate washout strength rhi = 
10-3 eV. In this case an intermediate behavior between those of the weak and 
strong washouts is observed: similarly to the strong washout case, the final 
asymmetry remains almost unchanged whether the CP scattering asymmetries 
are included or not. However, as in the case of weak washouts, its value remains 
well below what would be obtained by starting with an equilibrium density of 
N's. 



5 Conclusions 

In this paper we have computed the CP asymmetries of scattering processes in- 
volving the top quark and the gauge bosons. CP violation in scatterings gives 
an important contribution to the generation of a lepton asymmetry at high tem- 
peratures (T ^ 2Mi), and in particular in the zero temperature approximation 
adopted in our calculations this contribution is by far the dominant one. We 
have compared our results with the approximate expressions of the scattering 
CP asymmetries in terms of the decay asymmetry, concluding that in scenarios 
in which the heavy Majorana neutrino masses are sufficiently hierarchical this 
approximation provides reasonably accurate results. We have shown that when 
the sources of CP violation in scatterings are included in the BE, in the limit 
of very weak washouts a strong cancellation between the asymmetries generated 
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at early times (when < Y^^), and the asymmetry of opposite sign generated 
at later times (when Yj^ > Y^) takes place, sizably suppressing the final lepton 
asymmetry with respect to the cases in which CP asymmetries in scatterings are 
neglected. 

In the strong washout regimes, for which the final lepton asymmetry is almost 
independent of the conditions at early times, such as the initial value of the right- 
handed neutrino density, the final results are instead essentially unaffected by the 
inclusion of the new sources of CP violation. 
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6 Appendix: Three body decay 

The three body decay width for A^i — > iiQ^t is 

i ([(x - OQ + atf - ixat] [(1 - X - aef - 4xa^])'/' , (67) 

where ay = (m^/Mi)^ and ch ^ (F/z/Mi)^, with Th being the decay width of 
the Higgs boson. Fo(A^i iiH) = |AjipMi/167r is the two body decay width 
of A'^i (at zero temperature) and the integration variable is x = {pt +pq)'^/Mi, 
where pt and pg are the momenta of t and respectively. 

We assume rrit + > mn, which is generally valid at high T if thermal 
masses are considered and also at T = if the Higgs boson is not too heavy. 
In this case the Higgs boson exchanged in the internal line of the three body 
decay cannot be on-shell, and the Higgs width (parametrised by ch) can be 
neglected. (Note that a resonant contribution would in any case correspond to 
the two body decay A^i — > £iH rather than to a genuine three body process.) In 
the zero temperature limit at, aq, an Q the integral in eq. (f671) would get a 
large enhancement from the region corresponding to small values of x. However, 
for finite values of the thermal masses this enhancement is not present, and in 
particular for T/Mi > 10~^ the three body decay rate is always less than 6% of 
the two body decay rate. Note also that, due to the effects of thermal masses, 
the phase space for both decays actually gets closed when T approaches Mi [9] . 
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